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In this work, we show how the complete set of splitting functions relevant 
for the evolution of various distribution functions describing nucleonic helicity 
structure can be obtained in the light front Hamiltonian perturbation theory 
Q\ ■ using completely fixed light front gauge, A + = 0. 
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^j ! I. INTRODUCTION 

o ... 

How the total helicity of a nucleon is distributed among its constituents is an active area of 
research at present. In a previous work |j] Q we have shown that in a gauge fixed theory (with 
A + = gauge) total helicity operator can be satisfactorily defined and separated into orbital 
and intrinsic parts for quarks and gluons and defined the structure functions containing 
information regarding the orbital helicity distribution of the constituents. There we also 
employed our newly proposed methods of calculating dressed parton structure functions in 
light front Hamiltonian framework [TJ] and obtained some of the splitting functions which are 
relevant for the (Revolution of these structure functions and the corresponding anomalous 
dimensions, which were identical with the recent results [BJ using different techniques (but 
Ph| with the same A + = gauge). In this work, we shall use our gauge fixed formulation and 

find out the complete set of splitting functions relevant in the evolution of various helicity 
distributions to complete the picture. 

In our previous calculation of the various structure functions in paper I, we used the 
dressed quark and gluon targets which are eigenstates of P = (P + , P 1 , P 2 ) and S 3 . It was 
implicitly assumed there that the total P 1 = (P 1 , P 2 ) for the targets to be zero, which one 
could always do by appropriately orienting the 3-axis. This in turn forced the orbital motion 
of the dressed parton associated with its centre of mass motion to vanish. Consequently, 
our calculation produced the restricted number of splitting functions which are important 
for the QCD evolution of the assumed state. To obtain the complete set of relevant splitting 
functions and the corresponding anomalous dimensions || which, for example, takes into 
account how the orbital motion evolves under QCD interaction, one has to consider a target 
state where orbital motion of the constituents is present even in absence of QCD interaction. 
This is nothing but due to the centre of mass motion of the dressed parton. Note that this 
centre of mass motion has nothing to do with that of a composite system such as meson. For 
the calculation involving a meson target, it will be replaced by the internal motion of the 
concerned parton. We shall introduce the relevant target state for our calculation in Sec. II, 
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1 Henceforth referred, frequently in this work to avoid repetition, as paper I. 



where we have also introduced the multi-parton wave-functions relevant for our calculations. 
In Secs.III-IV we present our calculations and then conclude in Sec.V. 

II. THE TARGET STATE AND MULTI-PARTON WAVE-FUNCTIONS 

Required state for the dressed parton can be obtained by taking a suitable superposition 
of states where all the individual states carry definite longitudinal and transverse momenta 
(P) as well as the intrinsic helicity (a) and then projecting out those states which carry 
same j 3 by introducing a Kronecker-delta function. We consider the following as our target 
state 

I S) = EE e ™^ciL(^ + , I P 1 \,m)5 m>j 3-* \P,a), (1) 

a m 

where (f> parametrizes the state such that (P^P 2 ) = \P ± \(cos<j), sin</>). Note that we have 
written the ^-dependence separately as a sum over m, the eigenvalue of L 3 , so that $ cm is 
independent of 0. The dressed parton state \P,o) in the RHS of eq.(|l|) can be written in 
the Fock-basis using multi-parton wave- functions which are independent of momenta P [|J . 
For a dressed quark, it is given by 

dkfd 2 k±i dk 2 d 2 k± 2 



\P\ Pi, a) = y/J/Mk, a)\0) + ^f ^ d k± \ ^ k± \ ^nfP^5\P - kl - k 2 ) 



x $^ iA (P; h, k 2 )b\h, a 1 )a\k 2l A)|0) + •••}, , (2) 



where the normalization constant M q is determined perturbatively from the normalization 
condition, 

(P /+ , Pi, a'\P + , P ± ,a) = 2(27r) 3 P+^(P+ - P' + )S 2 (P ± - P|). . (3) 

Similarly for a dressed gluon 

dkfd 2 k±i dk 2 d 2 k± 2 



]P+ P ± , A) = yfiTJaHk, a)\0) + E / I \ T kX -\ ^(^fP^8\P -h-k 2 

V l £&J J2{2ixfkt J2(2n) 3 kt 



x^ 2 (P^i,k 2 p(k 1 ,(T 1 )d i (k 2 ,a 2 )\0) 

lr dkid 2 k±\ dk 2 d 2 k ±2 



+ _y i » ^ 2 ±2 h(27r) 3 P+S 3 (P -h- k 2 ) 

^tJ J2(2n) 3 ki J2(2n) 3 k+ y K ' K ' 



x $^ lA2 (P; h, k 2 )a\h, Xt)a\k 2 , A 2 )|0) + • • • } . (4) 



We introduce the boost invariant amplitudes ip° a 2 ( x > k± ) anc ^ so on by ^x 1 x 2 (P'i ^i> fa) 
-^ipa x 2 ( x > K± ) an d so on - Here the relative momenta (x, fir 1 ) are defined as 



k+ = xP + , k\ = k 1 + xP l and fc+ = (1 - x)P + , fcj = -k 4 + (1 - x)xP l . (5) 

From the light-front QCD Hamiltonian, to the lowest order in perturbation theory, we have 



rJ 2 9) (x,K ± ) = -S=T"^L^xt h^~ + -(o-i ■ K±)Ax*et , (6) 

/2(2tt) 3 «i I 1 - z x J 



<ep(*,*o = -^^f^.^j^, (7) 



^ K* A <M x "■>- "-L 1- X 



'"^i — ^ — t " K ^ + 7*° + r^ % / " Ai£a2£a • (8) 

Here we have already taken target mass and bare quark mass to be zero, M = m = 0, since 
nonzero masses in the above wave-functions are going to produce higher twist effects. Our 
next task is to calculate the structure functions with | S) being the dressed quark target 
[i.e., eq.(0) and eq.(|2|) combined] and dressed gluon target [i.e., eq.(|l|) and eq.(|4|) combined]. 

III. DRESSED QUARK STRUCTURE FUNCTIONS 

For the dressed quark eq. ([!]) can be written as 

\S) = e m ^$ cm (m 1 ) \P, + 1 -) + e"^$ cm (m 2 ) |P,~>, (9) 

where ni\ = j 3 — | and m 2 = j 3 + |. For such a state, in the zeroth order perturbation 
theory, i.e., if the two-particle wave-function in eq.@ is zero, Aq(x,Q 2 ) and Aqi(x,Q 2 ) 
come out to be proportional to 5-functions and are given by[| 

Aq(x, Q 2 ) = ^{ |$ cm (™i)| 2 - |$ cm (m 2 )| 2 } 5(1 -x) = ^AS 5(1 - x) , 

Aq L (x, Q 2 ) =M q {m l ^{m^ 2 + m 2 |$ cm (m 2 )| 2 } 5(1 - x) = L q 5(1 - x) , (10) 

while Ag(x,Q 2 ) and Agi(x,Q 2 ) are zero. Of course, Af q = 1 has been used above, since 
calculations are performed only in zeroth order in the coupling. 

Now, the structure functions can be calculated to the first order in light-front pertur- 
bation theory by truncating the dressed quark state in eq.(@) at the two particle level and 
using the wave- function in eq.(|6]). Straight-forward calculations give 

Aq(x, Q 2 ) = *k AS { 5(1 - x) + gc> In * i±^ } , (11) 

where A and /i are the upper and lower momentum cutoff. Note that it is same as we 
obtained in paper I except an extra weight factor AS coming from the fact that the target 
state now is a particular superposition. Similarly, 

a A 2 

Ag(l - x, Q 2 ) = N q AS -H-Cf In — (1 + x) . (12) 
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See paper I for the definitions of various structure functions. 



On the other hand, we get extra terms while calculating Aqi(x, Q 2 ) and Apx,(l — x, Q 2 ) as 
given below. 



<x« 



A 2 



Aq L (x, Q 2 ) = - A/" g AS -LCf In — (1 - x)(l + x 



2n 



f- 



KL q {5(l-x ) + ^C f l^^} 



(13) 



Ag L (l-x,Q 2 ) 



o: s 



A 2 



M„ AS — Cfln— x(l+ r 

A 2 



, g _ ^„. u2 



o. s 



+ N a L a ^C f \n— (l + o- 



•«~9 27r -^"„2 



(14) 



The extra terms in eq. (|C|) and eq. ( R ) are coming due to the fact that we have used |l. 



L q = -i 
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(15) 



where is defined by (k; 1 ,/^ 2 ) = |k Jl |(cos6', sin6>). In our previous calculation Jr did not 

contribute for we assumed P 1 - = 0. Also notice that the other terms in L q (as given in 
paper I) have no contribution for they produced odd integral of n 1 - over symmetric limit. 

Now, the normalization constant M q to the order a s can be calculated from eq.(^) and 
is given by M 



M,= {l-|i C/ l„^/ d *i±iL}, 

L 2ir u 2 J 1 - x > 



n ol s A 2 

J 2tc fi 2 



dx 



1 — x 
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2 



(16) 



Putting M q back into eqs. ([TT] - |T4D and keeping terms only to the order a s , we have 
Aq(x,Q 2 ) = iAS { 5(1 -x) + ^ln^ P S s( qq )(x) } , 
Ag(l-x,Q 2 )= AE^]n—P ss{gq) (l-x), 



2n /i 



a, , A 2 



Ag L (x,Q 2 ) = AS ^ln — P LS{gq) (x 

+ ^ { S(l 



a s , A 2 _ . , i 

- ,r) r — In — P L L{ qq )(x) }, 



2tt /x 2 



«, . A 2 



Ag L (l - x, Q 2 ) = AS -^ In — P LS{gq) (l - x 



2n fi 2 



LS( gq )\ 

a» . A 2 



+ L «^ ln ^ P LL{gq)\ 



l-x) 



(17) 



where we have defined the various splitting functions as follows. (Notice that the argument 
of glunic distribution functions is (1 — x) which is the momentum fraction carried by the 
gluon here.) 



PsS(qq){x) = Cf 


- 1 + x 2 


-6(1- 

2 v 


- x) 


PsS^i 1 ~ X ) = Cf(l + X) 






PLS(qq)(x) = ~Cf (1-X 2 ) 






P L S(gq){l ~ X) = -C f x(l + . 


*) 




PLL( qq )(x) = Cf 


- X 3 + X 

ki-x)+ + 


V 

2 v 


- x) 


PLL(gq){l-x) = 


Cf (1 + x 2 ) 







(18) 

Also, since there are no contributions coming from L q in Aq(x, Q 2 ) and Ag(x, Q 2 ), we have, 

PsL iqq) (x) = PsL igq) (x) = 0. (19) 

Last two splitting functions in eq. (|l8|) and that in eq. (|T9"|) are new compared to those in 
paper I. 

IV. DRESSED GLUON STRUCTURE FUNCTIONS 

Now we repeat the calculations given in the previous section by replacing the target state 
with a dressed gluon. For the dressed gluon eq.(|T]) can be written as 

\S) = e im ^ $ cm (mx) \P,+1) + e m ^ $ cm (m 2 ) \P,-1) , (20) 

where mi = j 3 — 1 and m 2 = j 3 + l- For such a state, in contrast to the quark case, Ag(x, Q 2 ) 
and Ag^(x,Q 2 ) are proportional to the 5-function in the zeroth order perturbation theory 
while Aq(x,Q 2 ) and Ag^(x,Q 2 ) are zero. 

Ag(x, Q 2 ) = N g { |$ cm (™i)| 2 " l$cm(m 2 )| 2 } 6(1 - x) = Ag 6(1 - x) , 

Ag L (x, Q 2 ) = Af g { rm |$ C m(mi)| 2 + m 2 |$ cm (m 2 )| 2 } 6(1 - x) = L g 6(1 - x) , (21) 

Now, with the dressed gluon state truncated at the two-particle level as given in eq.(|), 
we can calculate all the structure functions to first order in a s and the results are as follows 

a s A 2 
Aq(x, Q 2 ) = Ag -^-In— Pss( q9 )(x) 

j-i i\ Lb 

a s A 2 a s A 2 

Aq L (x,Q 2 ) = Ag -^ln— 2P LS(qg) (x) + L g -^ln— 2P LL(qg) (x) , 

A 2 

Ag(x, Q 2 ) = Ag { 6(1 - x) + ^ln— P SS (gg)(x) } , 



a s A 2 
Ag L (x, Q 2 ) = Ag -^-In— P LS ( gg )(x) 



2n ^ 

27T~'> 2 " 

6(-\ -x) + 

2tt /i 2 



A 9 

+ L g { 6(1 - X) + %Lln— PLL { gg)(x) } (22) 



where, like the quark case, we have defined the splitting functions as follows. 

PsS(qg)(x) = - [X 2 - (1-X) 2 }, 

, x (, ^ r l 1 , (1— X) 3 -| „, , /ll 1 



PLSto)(x) = 1 -[X 2 + {1-Xf}{l-X) 

PLS( ga )( x ) = 2N (x-l)(x 2 -x + 2), 

PLL(qg)(x) = - X[X 2 +(1-X) 2 ], 



2 

x 1 — X . . -| „,_ . /ll 



P LL{99) {x) = 2N x{ —^- + —±+x(l-x)} + 5(l- x) (-N - -) , (23) 

and 

P S L( qg) (x) = PsL ig9 )(x) = 0. (24) 

Here, in the above calculation, we have used the normalization constant J\f g as 



a. , A 2 



x 1 — x ,„ ,-, 1 



M g = 1 In — dx 2Nx{ 1 hi(l-i)} + - [ x 2 + (1 - x 

2ir [i l ' - -) 



1 — x x 2 

a s , A 2 r , 2N /ll T 1 

-±\n— dy- + (-N-- 

2n \i A J 1 — y v 6 3 



^ln^/dy^- + (^iV-J). (25) 



Last two splitting functions in eq. (p3|) and that in eq. (|24]) are new compared to those in 
paper I. Also, to the best of our knowledge, spin-dependent Altarelli-Parisi splitting func- 
tions denoted here as Pss(q g ) an d Pss( gg ) are explicitly derived using light-front Hamiltonian 
perturbation theory for the first time here. 

Corresponding anomalous dimensions are defined as 



A n 

A AB{ab) 



j dxx n ~ l P A B(ab)(x). (26) 



Using this definition the results given in Ref. || for the anomalous dimensions are exactly 
reproduced and it is needless to rewrite them here. Note that for these dressed parton states, 
expectation value of J 3 to this order in perturbation theory is independent of a s , as it should 
be due to kinematical nature of the helicity operator, 



2(2tt) 3 P+5 3 (0) 2 

(S\J 3 \S) 9 



2(2vr) 3 P+5 3 (0) 



Ag + L g . (27) 



Eqs.(|27|) explicitly verify the helicity sum rule for the two different states considered here 



and show the consistency of our calculation. 



V. CONCLUSION 

In this work, we have presented the calculations of the complete set of splitting functions 
relevant in the evolution of the various helicity distribution functions defined in the gauge 
A + = 0. This work in association with our previous work in paper I gives a satisfactory 
description of various parts of helicity distributions and their evolutions in the gauge fixed 
light front Hamiltonian formulation and in conformity with the work presented in Ref. ||. 

It should be noted that our formulation and results are obtained in a gauge fixed theory 
(same is true for the calculation in Ref. |J). At present, we have very little idea in which 
processes the orbital helicity distributions that we have defined can be measured experimen- 
tally. It needs to be admitted that in absence of the experimental support and/or our results 
are reproduced by some other means, the gauge fixed formulation is little unsettled. On the 
other hand, gauge invariance in separating the helicity operator into intrinsic and orbital 
parts for quarks and gluons is a long standing issue. Lot of work has been done in the last 
few years towards defining such operators showing full respect to the gauge invariance [Q. 
Also, there are attempts to define operators compatible to the residual gauge invariance in 
light-front gauge ||. In spite of all these works, the issue still seems to be a little unsettled 
|7]] as well. In view of this dilemma, it needs thorough investigation to settle down the ques- 
tion of gauge invariance and/or find out processes sensitive to OAM that can be measured 
in the experiment. We hope to undertake such investigations in future. 
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APPENDIX A: 

In this appendix, we provide the necessary details of calculation in arriving at the ex- 
pression for AqL(x,Q 2 ) given eq.fll3|) starting from its definition as an example. All the 
other splitting functions can be obtained in a similar manner. We have (see paper I), 

&q L (x, Q 2 ) = — ^ J d V e~ irix (S | ^{Ol+i^d 2 - x 2 d 1 )^(0) + h.c] | S) (Al) 

with i] = iP + £~. Note that the operator involved here is a number conserving operator, 
so that the non-vanishing contributions will only be diagonal in particle number. Let us 
first consider the contribution coming from the first term in \S) given in eq.(^). Using 
the Fourier expansion for the dynamical fermionic field ifj + and the standard commutation 
relation between creation and annihilation operators (see paper I), we get two non- vanishing 
contributions. First one coming from the one particle sector and is simply given by 



Aq L (l)= M q l^cmK)! 2 e~ im ^(-i|-)e imi * 5(1 



x) 



M q m 1 \^ cm (m 1 )\ 2 5(l-x). (A2) 



The second one coming from the two particle sector is given by 
Ag L (2)= X q \® cm (mi)\ 2 £ /c/V e" imi * < A /(x, 

im-i_(j> 



K 
o"i,A2 
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= Ag$ cm (™i)l 2 E Jd 2 ^Ui-^l lX ;(x,^)(-z^\^ lX2 (x,^) 

+ A^m 1 |$ cm (m 1 )| 2 J2 / rf V ^^l^i^^W^x^,^) (A3) 

o"i,A2 

Here and in the following f implies helicity +| for quarks and +1 for gluons and so on. 
Using xi — (1)0), X-i — (0,1) and e±i = -4r(l,±i) in eq.(H), various components of the 
two-particle wave-function can be simplified and are given by 



2(27r) 3 l«- L IVl-a; ^2(2^) 3 l^l\/l-^ 

^ T (x,« J -) = Vl i (x,« J -) = 0. (A4) 

Last two components are zero since we have assumed m q = 0, which forces the quark helicity 
flip interaction to vanish. Putting them back into eq. ( |A3| ) and performing k -integration, 
we get the required result from the two-particle sector. Thus the total contribution from 
the first term of \S) becomes 

A 2 

Aqi(l) + Aq{(2) = -Ar q \<S> cm (mi)\ 2 ^C f \n-(l-x)(l + x) 

Z7T fj, 

+ Af q m 1 \^ cm (m l )\ 2 {S(l-x) + -^C f \n—\^}. (A5) 

— i\ ft i- Ju 



The contribution from the second term in \S) is again given by eq.(|A2|) and eq. (|A3|) with 
mi —>■ m,2 and |— >|. Then one uses 



WllfaK ) = - r-— 1^1 /I = ' ni( X ' K ) = - r-^,. n /q = » ( A6 ) 



ip^ix,^) = i) l u (x,^) = 

gT a V2xe- ie , ^ _ ^ yg, 

^2(2^)3 I^IVl- a:' U X ' " ^2(2tt) 3 I^IVl - ^ 

to obtain 

r> A 2 

Aq I L I (l) + Aq{ I (2) = Af q \$ cm (m 2 )\2 -^ Cf ]n— (l-x)(l + x) 

r rv A t -4- t i 

+ ^ m 2 |$ cm (m 2 )| 2 { 5(1 -^H-^C/ In— ^^ } . (A7) 

Note that eq. (|A~5|) and eq. (^71) are similar except the first term having opposite sign and m\ 
is replaced by m%. Adding them together, we obtain the final result as given in eq . (p~3|) . 
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